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Abstract 



o 

■ For the analysis of systems consisting of small, regular objects, the methods of math- 

^ . ematical morphology applied to images of these systems are well-suited. One of these 

methods is the use of Voronoi polygons. It was found that the Voronoi tessellation method 
represents a powerful tool for the analysis of thin film morphology and provides nanostruc- 
tural information to many multi-particle assemblies. In these notes, several morphological 
\^ [ algorithms are analyzed and we study how to join all of them to design a graphical user 

Tij" ' interface (GUI) that provides as input for the system the "AFM image" and interprets the 

^] ■ output of the system in terms of errors and generators coordinates. 
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X ■ 1 Introduction 



In the last decades, there has been an increasing interest in a geometrical construct called the 
Voronoi diagram (e.g. [1], [2], [3] and [4]). The Voronoi diagram is a data structure extensively 
investigated in the domain of computational geometry (e.g. [5]). 

Given some number of points in the plane, their Voronoi diagram divides the plane according 
to the nearest-neighbor rule: Each point is associated with the region of the plane closest to it, 
so it is a tessellation of M^. We have already noted that the concept of the Voronoi diagram 
is used extensively in a variety of disciplines and has independent roots in many of them (e.g. 
[6]). The first extension of them was to the area of crystallography (the area we are interesting 
in), works in this field are for example [7] and [8]. 

Since there is a large number of empirical structures which also involve tessellations of M^, 
one of the most direct applications of Voronoi concepts is in the modelling of such structures 
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and the processes that generate them. In these notes, we use the Voronoi assignment model 
in the modelhng of physical-chemical systems. Such systems under study consist of a set of 
sites occupied by atoms, ions, molecules, etc. (depending on the specific application) which are 
represented as equal-size spheres. Our system is formed by sites regularly arranged in M^, they 
assume form of lattice (the structure is said to be crystalline). Thin metal films images with 
Atomic Force Microscopy (AFM) consist of small two-dimensional islands (objects) distributed 
on the substrate. The quantitative characterization of the object arrangement can bring in- 
formation about internal processes in the studied system. We apply methods of mathematical 
morphology to thin metal films images with Atomic Force Microscopy, to assign the model: 
The Voronoi Growth Model. Voronoi polygons has been employed for providing nanostructural 
information to these multi-particle assemblies. We analyze morphological algorithms applied 
to these tessellations, e.g. to restore the generators from a given Voronoi diagram. 

As a graphical user interface (GUI) makes easier for the user to obtain information from 
algorithms, we present how to join all algorithms, we have studied, to design one. The graphical 
user interface provides as input for the system the "AFM image" , and interprets the output in 
terms we are interesting in. We note that this work can easily be extended, if we have images 
from other fields like ecology, meteorology, epidemiology, linguistics, economics, archeology or 
astronomy, that we suspect are a Voronoi diagram. 

The structure of this paper is as follows. In Section 2 we recall the mathematical theoretical 
background about Voronoi diagrams and we give an application of them, that it is called the 
"Voronoi Growth Models "which we will use in the analysis of film nanographs. In Section 3 
we explain the mathematical solution to the problem proposed here. Next in Section 4, we give 
and analyze algorithms of the mathematical solution and in Section 5 we finish with important 
concluding remarks and directions for further research. 

2 Preliminaries 

In this section we will review the basic notions we shall require for the sections to follow. For 
more details about them we refer to [9] and [10] for the first investigation of mathematical as- 
pects of Voronoi diagrams, [11] and [12] for papers that present surveys about Voronoi diagrams 
and related topics, and [13] for a good introduction to all applications of Voronoi diagrams to 
sciences. 

2.1 Mathematical Background 

We will define the Voronoi diagram and introduce properties and notations to be commonly 
used in this notes. 

We work with a finite number, n, of points in the Euclidean plane, and assume that 
2 < n < oo. The n points are labeled by pi, ■ ■ ■ , p„ with the Cartesian coordinates (xn, X12), ■ ■ ■ , 
(a;„,i,x„2) or location vectors xi, ■ ■ ■ ,x„. The n points are distinct in the sense that Xj 7^ Xj 
for i 7^ j, i,i e /„ = {1, ■ ■ ■ Let p be an arbitrary point in the Euclidean plane with 
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coordinates {xi,X2) or location vector x. Then the Euchdean distance from p to Pi is given by 

d{p,Pi) =11 X - Xi 11= a/ (Xi - XiiY + (X2 - Xi2y. 

If Pi is the nearest point from p or is one of the nearest points from p, we have the relation 
II X — Xj ||<|| X — Xj II for j ^ i, j ^ In- In this case, p is assigned to pi. Therefore, 

Definition 2.1. Let P = {pi, ■ ■ ■ ,Pn} where 2 < n < oo and Xj 7^ Xj for i 7^ j, i,j e /„• We 
call the region given by 

y(pi) = {x s.t. II X - Xi ||<|| X - Xj II for j 7^ 2, j e /„} (1) 

the (ordinary) Voronoi polygon associated with pi (or the Voronoi polygon of pj), and the set 
given by 

V = {V{pr)r-- .V{pn)} 
the (planar ordinary) Voronoi diagram generated by P (or Voronoi diagram of P). 

We can extend the above definition to the m-dimensional Euclidean space, but for our 
proposes we only need the Euclidean plane. So, we shall often refer to a planar ordinary 
Voronoi diagram simply as a Voronoi diagram and an ordinary Voronoi polygon as a Voronoi 
polygon. 

For a Voronoi diagram V we have the following definitions. 

Definition 2.2. We call the Pi of V{pi) the generator point or generator of the ith Voronoi 
polygon, and the set P = {pi, ■ ■ ■ ,p„} the generator set of the Voronoi diagram V (Figure 1). 

Notation 2.3. For brevity we may write for V{pi). Also we may use V{xii,Xi2) or V^(xj) 
when we want to emphasize the coordinates or location vector of the generator pi. In addition, 
we may use V(P) when we want to explicitly indicate the generator set P of V. 




Figure 1: Voronoi diagram 
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Definition-Notation 2.4. Given a Voronoi diagram V, since a Voronoi polygon is a closed 
set, it contains its boundary denoted by dV{pi). The boundary of a Voronoi polygon may 
consist of line segments, half lines or infinite lines, which we call Voronoi edges. Noticing that 
= is included in the relation of equation ([1]), we may alternatively define a Voronoi edge as a 
line segment, a half line or an infinite line shared by two Voronoi polygons with its end points. 
Mathematically, if V{pi) fl Vpj) ^ 0, the set V{j)i) fl Vpj) gives a Voronoi edge (which may be 
degenerate into a point). We use e{pi,pj) for V{pi) fl Vpj), which is read as the Voronoi edge 
generated by Pi and pj. Note that e{pi,pj) may be empty. If e{pi,pj) is neither empty nor a 
point, we say that the Voronoi polygons V{pi) and V{pj) are adjacent. 

An end point of a Voronoi edge is called a Voronoi vertex. Alternatively, a Voronoi vertex 
may be defined as a point shared by three or more Voronoi polygons. We denote a Voronoi 
vertex by qi (see Figure 1). When there exits at least one Voronoi vertex at which four or 
more Voronoi edges meet in the Voronoi diagram V, we say that V is degenerate (Figure 2); 
otherwise, we say that V is non-degenerate. 



In the previous definitions of Voronoi diagram V, we have defined a Voronoi diagram in an 
unbounded plane. In practical applications, however, we often deal with a bounded region S, 
where generators are placed. In this case we consider the set given by 



We observed that an ordinary Voronoi diagram consists of polygons, as a polygon can be 
defined in terms of half planes, we have the equality of Proposition 12.61 

Notation 2.5. Given a Voronoi diagram V(pi, ■ ■ ■ ,Pn), we consider the line perpendicularly 
bisecting the line segment pip] joining two generators pi and pj. We call this line the bisector 
between pi and pj and denote it by b{pi,pj). 

Since a point on the bisector b{pi, pj) is equally distant from the generators Pi and pj, b{pi,Pj) 
is written as 



*■ 



* 



Figure 2: Degenerate Voronoi diagram 



vns = {V{pi) n 5, ■ ■ ■ , 



v{Pn) n S} . 



KPi^Pj) = S.t. II X - Xj 



x-Xj ID, j ^ i. 
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The bisector divides the plane into two half planes and gives 

H {pi, pj) = {:si s.t. II X - Xi ||<|| X - 11} , jj^i. 
We call H{pi,pj) the dominance region of pi over pj. 

Proposition 2.6. Let P = {pi, ■ ■ ■ C M^, where 2 < n < oo and Xi ^ xj for i ^ j , and 

i,j G In- Then 

Vip.)= fl Hip,,p,) 
ie/n\{i} 

where V{pi) is the (ordinary) Voronoi polygon associated with pi and set 

ViP) = {V{p^),--- ,V{pn)} 

where V{P) is the (planar ordinary) Voronoi diagram generated by P. 

As a degenerate Voronoi diagram requires special lengthy treatments which are not essential 
we avoid this difficulty and we often make the following assumption: 

Assumption 2.1.1. (the non- degeneracy assumption) . Every Voronoi vertex in a Voronoi 
diagram is incident to exactly three Voronoi edges. 

The largest empty circle in a Voronoi diagram. 

Definition 2.7. For a given set P of points, if a circle does not contain any points of P in its 
interior, the circle is called an empty circle. 

Theorem 2.8. Let Q = {qi,--- ,qnv} be the set of Voronoi vertices of a Voronoi diagram 
generated by P. For every Voronoi vertex, Qi & Q, there exists a unique empty circle Ci centered 
at qi which passes through three or more generators. Under the non- degeneracy assumption, Ci 
passes through exactly three generators (Figure 3). 

From this theorem, the non-degeneracy assumption (Assumption 12.1.11) is equivalent to the 
following assumption. 

Assumption 2.1.2. (the non co-circularity assumption) Given a set of points 
P = {pi, ■ ■ ■ ,pn} C (4 < n < oo), there does not exist a circle, C , such that pn, - ■ ■ ,Pik, 
k > 4, are on C, and all points in P\ {pn, ■ ■ ■ ,Pik} are outside C. 

Theorem 2.9. Circle Ci in Theorem \2.8{ is the largest empty circle among empty circles cen- 
tered at the Voronoi vertex qi . 

Here we have seen some properties of the Voronoi diagram, but it has many more. For 
example, if one connects all the pairs of sites whose Voronoi cells are adjacent then the resulting 
set of segments forms a triangulation of the point set, called the Delaunay triangulation. 
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Figure 3: Empty circle centered at a vertex of a Voronoi diagram 
2.2 Voronoi Growth Models 

Thin metal films deposited on a surface consist in their initial stage of growth of small islands. 
Basic information about nucleation processes during the thin film growth can be derived by the 
morphological analysis of the film AFM image. For very thin metal film or generally for systems 
consisting of small regular objects, the methods of mathematical morphology are well-suited 
to the study of spatial distribution of objects in images (e.g. [14]). There is a large number 
of empirical structures which involves tessellations of (and more generally in M™), one of 
the most obvious direct applications of Voronoi concepts is in the modelling of such structures 
and the processes that generate them. These models produce spatial patterns as the result of 
a simple growth process with respect to a set of n points (nucleation sites), P = {pi, ■ ■ ■ ,Pn}, 
at positions xi, ■ ■ ■ , x^, respectively, in or a bounded region of M^. If we make the following 
assumptions, the resulting pattern will be equivalent to the ordinary Voronoi diagram V(P) of 
P: 

Assumption 2.2.1. Each point pi (i = 1, ■ ■ ■ ,n) is located simultaneously. 

Assumption 2.2.2. Each point Pi remains fixed at Xi throughout the growth process. 

Assumption 2.2.3. Once pi is established, growth commences immediately and at the same 
rate U in all directions from pi. 

Assumption 2.2.4. /j is the same for all members of P. 

Assumption 2.2.5. Growth ceases whenever and wherever the region growing from pi comes 
into contact with that growing from pj (j ^ i). 

Together, Assumptions I2.2TII2.2.5I define the Voronoi Growth Model. The Figure 4 shows a 
series in stages in such a growth process. 

Generally speaking, one obvious application is to model crystal growth about a set of nu- 
cleation sites. Here assumptions l2.2.1tlXT^ are equivalent to assuming an omni directional, 
uniform supply of crystallizing material to all faces of the grind crystal in the absence of any 
absorbable impurities. Assumption 12.2.31 also implies that the rate of growth of the volume of 
a crystal will be proportional to its surface area. 
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Figure 4: Growth process of crystals (Courtesy of Pablo Stoliar) 



Also, growth models to modelling phase transitions in metallurgy involving transformation 
of an isotropic, one-component solid through nucleation, and isotropic growth of grains of a 
new or re-crystallized phase. In this context the Voronoi growth model is sometimes referred 
to as the cell model or the site saturation model. Specific examples include the covering of a 
metallic surface by films or layers of corrosion product where the nucleation sites P might be 
surface imperfections such as impurities, points of intersection with bulk defects and surface 
pits. Another example is the growth of thin films of metal or semiconductors. In these examples 
if the thickness of the film is small relative to the spacing between the nucleation sites or if the 
grain boundaries are perpendicular to the plain of the film, a two-dimensional representation 
is appropriate. 

3 Problem-Solution 

If the Voronoi assignment and growth models described in the previous section are appropriate 
for modelling a particular phenomenon, we would expect spacial patterns of the phenomenon 
to display characteristics of Voronoi diagrams. In case we have a tessellation S, we have to 
consider ways of determining if it is a Voronoi diagram based on some set P (this problem 
has been studied e.g. in [15]). Recognizing a Voronoi diagram is closely related with the next 
generator recognition problem. 

The generator recognition problem: Provided that the Voronoi edges of a non-degenerate 
Voronoi diagram V(P) are given, we recover the locations of generators P. 

3.1 The generator recognition problem 

The first problem we approach is to restore the generators from a given Voronoi diagram, that 
is, the inverse problem of constructing the Voronoi diagram from the given points. For this 
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problem itself, we propose the following geometrical approach (see e.g. [13]). 



Let Qi be a Voronoi vertex, Pii,Pi2,Pi3 be generators whose Voronoi polygons share g^, and 
Qii, (li2, Qis be the Voronoi vertices of the Voronoi edges incident to g,. From Theorem 12.81 qi is 
the center of the circle that passes through Pii,Pi2,Pi3- Since Voronoi edges e(pii,Pi2), ^{j>i2,Pi3) 
and eipa.Pii) perpendicularly bisect line segments Pi2Pi2, Pi2Pi3, PisPn, respectively, we have 
the equations: 



Hence 2ai + 2/3j + = 27r, i.e. = vr — /3j — 7^ = vr — Zqi2qiqi3. From this equation, we 
obtain the following theorem. 

Theorem 3.1. Let qi/jj be a Voronoi edge in a non- degenerate Voronoi diagram, and 6jk, 
k = 1,2,3 be the acute angles at qi and qj, respectively, where k is indexed counterclockwise 
from qiqj at q^ and clockwise at qj. Let Lik (Ljk) be the half line radiating from q^ (qj) with 
angle ir — 9i2 fvr — 9j2) with qiqj in the sector of Oik ((^jk), k = 1,3 . Then the intersection 
point made by Ln and Lji, and that by Li^ and Lj^ give the generators of the Voronoi diagram 
sharing q[qj. 

We develop this theorem into a more general theorem with which we can examine whether 
or not a given planar tessellation S = {Si, ■ ■ ■ , Sn}, is a Voronoi diagram. We suppose that 
the tessellation S consists of convex polygons and every vertex has exactly three edges. Let 
Qiir ' ' ) Qiki be the vertices of a polygon Si indexed counterclockwise. Let pij be the intersection 
point in Si obtained through the same procedure stated in the previous theorem, where g^gj is 
replaced by qiJqiP^, j = 1, ■ ■ ■ ,ki {ki + 1 should be read as 1). Then we have the main theorem. 

Theorem 3.2. A planar tessellation consisting of convex polygons whose vertices are all degree 
three is a Voronoi diagram if and only if Pn = Pi2 = ■ ■ ■ = Piki holds for i G /„„ , where pij is 
defined in the above. 

An alternative method to find the generators P = {pi, ■ ■ ■ ,pn} of a Voronoi diagram V 
proposed in [16] is the following algebraically method. It is based on the perpendicular bisector 
property i.e. the line segment pkpi joining the generators pk and pi of two adjacent Voronoi 
polygons V{pk) and V{pi) is bisected perpendicularly by the common edge e{pk,pi) of V{pk) 
and V{pi) (see Proposition 12. 6p . So, that means that p^ and pi are subject to the following 
conditions: 

Condition 1. p^ and pi lie on a line perpendicular to e{pk,Pi). 

Condition 2. pk and pi are equidistant from e{pk,pi). 

These conditions can be formulated algebraically to form a linear system of equations which 
can be solved to find the locations pk and pi. Let (si,S2) and (ti,t2) be the locations of the 



Zpiiqiqn 
Z.piiqiqi2 
^Pi2qiqi3 



Z.pi-iqiqii = ai, 

Zpi2qiqi2 = A, 
^Pi3qiqi3 = li- 
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end points of the common edge e{pk,Pi) of two adjacent members Vk and of V. We search 
{xk, Vk) and (x;, yi) the locations of the generators pk and p/, respectively. The segment e{pk,pi) 
lies on the line 



(S2 - t2)X - (Si - h)Y + t2(Sl - tl) - h{s2 -t2)=0 

As V is a Voronoi diagram, Condition 1 gives 



(xfc - xi){si - ti) + {yk - yi){s2 - ^2) = 



and Condition 2 gives 



■52 ~ ^2 

Si -ti 



)xk -yk + { 



S2 — ^2 
Si - ti 



)xi - yi 




S2 ~ ^2 

Si -ti 




Suppose that V is non-degenerate and has m interior edges. Condition 1 gives a system of 
m equations and 2n unknows, and Condition 2 gives another system of m equations and 2n 
unknows. Taken jointly all equations we have enough constraints to provide a least squares 
solution for {xi,yi,X2,y2, ■ ■ ■ ,Xn,yn)- Specific methods for solving this equations are given in 
[16]. The point is that if V is a Voronoi diagram, then all equations will yield the same solution 
for (xi, yi, X2, y2, - ■ ■ , Xn, yn)- Evans and Jones in [16] outline three algorithms for its solution. 
Unfortunately, the algorithms require the inversion of poorly-conditioned matrices and may 
thus be highly unstable. 

3.2 Approximating Voronoi diagrams 

If the exact Voronoi diagram were given, we could determine the position of the generators by 
the previous methods of Subsection 13.11 However, such a situation is unrealistic. Recognizing 
that the recording of many types of empirical patterns often involve some measurement error, it 
is usual that given a pattern could not correspond to a Voronoi diagram even when we suspect 
that the pattern was generated by processes such as those in the Voronoi diagram. Even if 
theoretical consideration tells us that the diagram which appears in a phenomenon should be a 
Voronoi diagram, the error in observation process must perturb the original diagram. Therefore, 
the geometrical method would always tell us that the diagram is not the Voronoi, i.e., it would 
give us no information in almost every case. Methods proposed in [16], [17], [18] tell us at least 
approximate positions of the generators. 

4 Voronoi diagram in Computational geometry: Algo- 
rithms 

Geometric objects such as points, lines, and polygons are the basis of a broad variety of im- 
portant applications and give rise to an interesting set of problems and algorithms. Computers 
are being used to solve larger-scale geometric problems. Computational Geometry has been 
developed as a set of tools and techniques that takes advantage of the structure provided by 
geometry. 
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Now, we describe algorithms to solve the generator recognition problem (see [PFL]). First 
of all we must store the tessellation of which we are seeking the generators. A tessellation is 
typically stored as a list of vertex coordinates and its associated contiguity lists: lists which 
provide, for each vertex, the indices of the other vertices to which it is connected. If a vertex 
lies on an infinite edge, we store both the vertex and an arbitrary other point p on the infinite 
edge, where p is labeled a dummy vertex and given no adjacency list. In the input of the 
algorithms we will describe, we require that the number of ordinary vertices and the number of 
dummy vertices be specified, and that the dummy vertices be placed at the end of the vertex 
list (we will consider them as degenerate). 

Let V be a tessellation of the Euclidean plane and Q = {qi, ■ ■ ■ ,g„„} the set of vertices 
in which the last Uc vertices lies in a infinite edge. By Theorem 13.11 if the tessellation is a 
Voronoi diagram for each vertex qj of a given polygon Vi we can define a half line Lij of a giving 
direction radiating from qj into the interior of V^. The intersection of any two such half lines 
gives the location of the generator of Vi. 

[19] gives the implementation of the following algorithms. 

Algorithm I (naive) 

The above introduction suggests the following naive algorithm for tessellations such that all of 
whose polygons contain at least two non degenerate vertices: 

ALGORITHM I 

Step 1. Specify the polygons Vi, - ■ ■ ,Vn. 
Step 2. For each polygon V^: 

2.1 Find any two non degenerate vertices outlining Vi, say qii,qi2- 

2.2 For each vertex qik {k = 1,2) find the ray extending from qik {k = 1,2) through the 
generator in Vi, as we described above. 

2.3 Find the intersection of this two rays. 
PROBLEMS 

i) The requirement that each cell contain at least two non degenerate vertices. 

ii) Only two rays are used to determine the generator in each polygon. 

iii) If the two rays in a polygon are perfectly parallel (a simple modification is to find an 

additional ray emanating from a different non degenerate vertex in the polygon). 
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Algorithm II 

Errors in the generator determination of the previous algorithm could be minimized by using 
all the available rays rather than just two. Hence an alternative is the following algorithm: 

ALGORITHM II 

Step 1. Specify the polygons Vi, ■ ■ ■ , Vn. 
Step 2. For each polygon Vf. 

2.1 Find all non degenerate vertices outlining Vi. 

2.2 For each vertex find the ray associated. 

2.3 Find the intersection of every possible pair of rays. 

2.4 Average these intersection points. 

PROBLEMS 

i) The generator location errors using the Algorithm II are in fact typically considerably larger 
than for Algorithm I!!! 

Algorithm III 

This increase in error in the previous algorithm is attributable to the instability in intersecting 
certain select pairs of rays, one may modify Step 2.4 of Algorithm II by computing a weighted 
average of the intersection points, weighting each point according to an estimate of its stability, 
as in the following algorithm: 

ALGORITHM HI 

Step 1. Specify the polygons Vi, ■ ■ ■ , Ki. 
Step 2. For each polygon V^: 

2.1 Find all non degenerate vertices outlining V^. 

2.2 For each vertex find the ray associated. 

2.3 Find the intersection pik of every possible pair of rays in the polygon Vi. 

2.4 For each pair [k, I) of rays in the polygon Vi, estimate the stability of its intersection 
by perturbing the slopes of each of the rays by a small amount in either direction 
and seeing how much the intersection point changes. Record 6k^i = the sum of the 
sizes of these changes. 

2.5 Gompute a weighted average of the intersection points, giving pi^. the weight 

(4,/)-V(Sfc',/'(4',/')"')- 
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Note that a potential alternative to algorithms II and III is to find the point minimizing 
some penalty function such as the sum of squared perpendicular distances to the rays. The 
(weighted) averaging in algorithms II and III is equivalent to finding the location minimizing 
the (weighted) sum of squared distances to the intersection points of the rays. 

Remarks 

Algorithms I, II, and III are all entirely local; each polygon is determined solely based on its 
own vertices and their neighboring vertices. The accuracy of the algorithms can potentially be 
improved by incorporating information from neighboring polygons, e.g. by using the perpendic- 
ular bisector relation of Proposition 12.61 Paik, Ferguson and Li suggest modifying Algorithms 
I, II, and III to improve the results. 

All of the algorithms proposed are extremely fast, requiring just 0{n) observations, where n 
represents the number of generators to be determined. 

The errors in the inversion algorithms proposed are very small. However, in [19] they inquire 
about the size of the errors resulting when one of the vertices is recorded substantially in error. 



5 Summary and proposed extensions 



We start with an AFM image, like one in the figures below that represent the growth of crystals 
with different velocities (Courtesy of Pablo Stoliar). 
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They can be represented as a tessellation of a bounded region of the plane. We want to 
apply Voronoi diagrams to pattern recognition to this branch of solid state physics. If we 
suspect that these structures are generated by spatial processes resulting in tessellations which 
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can be constituted by Voronoi diagrams, we would have to follow the next steps to analyze the 
images in order to obtain properties of the thin films they represent. 

1. We approximate the image to extract the vertices of the tessellation. 

2. Using one of the algorithms of Subsection 3.2 we approximate this tessellation to a Voronoi 
diagram. 

3. In order to apply algorithms of Section 4, we store the vertices like we said in this section. 
The algorithms take as input this list of vertices and their adjacency lists, the output give 
us the generator points of the Voronoi diagram. 

4. The last step will be to measure the errors (root-mean-squared errors in the vertices 
locations) between the tessellation we obtain directly from the image, and the Voronoi 
diagram we obtain with the outputs of the algorithms I, II and III. 

Finally, the present work may be extended in the following way: every step uses a different 
computational algorithm, it is interesting to join all this steps to design a graphical user interface 
(GUI) that provides as input for the system the "AFM image" , and interprets the output of the 
system in terms of errors and generators coordinates. A user interface makes easier for the user 
to interact with the designed programs utilizing toolbar buttons and/or icons. Every software 
package, the one for step 1, step 2, step 3 and step 4, need a graphical user interface design that 
can be developed, we think it would be useful to do only one graphical user interface with all 
these algorithms inside, in that way we can obtain quickly the thin film information we need. 

This graphical user interface could be easily extended to another kind of images, i.e. we can 
provide as input for the system a variety of images that it could be represented by a Voronoi 
diagram. So, this work will be valuable not only in the field of crystallography, but also in 
the fields such as ecology, meteorology, epidemiology, linguistics, economics, archeology and 
astronomy where Voronoi diagrams are applied. 
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